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We discuss a dynamic procedure that makes the fractional derivatives emerge in the time asymp- 
totic Umit of non-Poisson processes. We find that two-state fluctuations, with an inverse power-law 
distribution of waiting times, finite first moment and divergent second moment, namely with the 
power index jj, in the interval 2 < fi < 3, yields a generalized master equation equivalent to the 
sum of an ordinary Markov contribution and of a fractional derivative term. We show that the 
order of the fractional derivative depends on the age of the process under study. If the system is 
infinitely old, the order of the fractional derivative, ord, is given by ord — 3 — fi. A brand new 
system is characterized by the degree ord = fi ~ 2. If the system is prepared at time — ta < and 
the observation begins at time t = 0, we derive the following scenario. For times < t << ta the 
system is satisfactorily described by the fractional derivative with ord — 3 — n. Upon time increase 
the system undergoes a rejuvenation process that in the time limit t » ta yields ord — fi — 2. The 
intermediate time regime is probably incompatible with a picture based on fractional derivatives, 
or, at least, with a mono-order fractional derivative. 

PACS numbers: 05.40.Fb, 05.60.Cd, 02.50.Ey 



I. INTRODUCTION 

The fractional calculus has recently received a great 
deal of attention in the physics literature, through the 
publications of books review articles |1,|3|, as well 

as an ever increasing number of research papers, some of 
which are quoted here The blos- 

soming interest in the fractional calculus is due, in part, 
to the fact it provides one of the dynamical foundations 
for fractal stochastic processes nil. The adoption of the 
fractional calculus by the physics community was inhib- 
ited historically because there was no clear experimental 
evidence for its need. The disciplines of statistical physics 
and thermodynamics were thought to be sufficient for de- 
scribing complex physical phenomena solely with the use 
and modifications of analytic functions. This view was 
supported by the successes of such physicists as Lars On- 
sager, who through the use of simple physical arguments 
was able to relate apparently independent transport pro- 
cesses to one another, even though these processes are 
associated with quite different physical phenomena [l3l| . 
His general arguments rested on three assumptions: 1) 
microscopic dynamics have time-reversal symmetry; 2) 
fluctuations of the heat bath decay at the same rate as 
do macroscopic deviations from equilibrium and 3) phys- 
ical systems are aged. We refer to assumption 2 as the 
Onsager Principle and show that it is tied up with as- 
sumption 3. 

Onsager's arguments focused on a system that is in 
contact with a heat bath sufficiently long that the bath 
has come to thermal equilibrium and consequently the 
system is aged. In statistical physics we know that the 
bath is responsible for both fluctuations and dissipation, 
and if the fluctuations are white the regression of per- 



turbations of the bath to equilibrium is instantaneous. 
This means that the energy absorbed from the system 
of interest by the bath, through macroscopic dissipa- 
tion, is distributed over the bath degrees of freedom on 
a very much shorter time scale than the relaxation time 
of the system. This property is summarized in the well 
known fluctuation-dissipation theorem, which has even 
been generalized to the case where the fluctuations in 
the bath do not regress instantaneously 0|. 

The dynamics of the physical variables to which the 
Onsager Principle apply are described by two different 
kinds of equations: 1) the Langevin equation, a stochas- 
tic differential equation for the dynamical variable and 
2) the phase space equation for the probability density. 
Two distinct methods have been developed to describe 
the phase space evolution of the probability density: the 
master equation introduced by Pauli and the continuous 
time random walk (CTRW) approach of MontroU and 
Weiss 15]. The CTRW formalism describes a random 
walk in which the walker pauses after each jump for a 
sojourn specified by a waiting time distribution function. 
It was shown by Bedeaux et al. that the Markov 

master equation is equivalent to CTRW if the waiting 
time distribution is Poissonian. However, when the wait- 
ing time distribution is not exponential, the case we con- 
sider here, the equivalence between the two approaches is 
only maintained by generalizing to the non-Markov mas- 
ter equation, the so-called GME Recently, Metzler 
argued that the GME unifies the fractional calculus 
and CTRW. 

Allegrini et al. have shown that to create a master 
equation compatible with the Onsager Principle requires 
that the system be entangled with the heat bath in such 
a way that the bath does not regress to equilibrium in- 
finitely fast. The system-bath entanglement is the result 
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of a rearrangement process that may take infinitely long 
to complete, leading to t he g eneralized master equation 
(GME) of Kenkre et al. [l2|. Allegrini et al. were 
concerned with how to make the GME stationary and 
therefore compatible with the Onsager Principle. Herein, 
we extend that discussion to include the connection with 
both the fractional calculus and the nonstationary con- 
dition. 

For simplicity we restrict our discussion to the case 
of a two-site system coupled to a heat bath, a problem 
closely related to the quantum problem of decoherence 
of a qubit, due to coupling to the environment. The con- 
nection of the CTRW, a two-state non-Markov master 
equation and the fractional calculus has been explored 
by Sokolov and Metzler Herein we go beyond their 
discussion and investigate the relationships among the 
fractional calculus, a particular generalization of the On- 
sager Principle, and one concept of aging. We also show 
that the result of Sokolov and Metzler, which refers to 
the young state, is an attractor for all the systems that 
are partially aged, and not infinitely aged. We call this 
process rejuvenation. 



A. Beyond the Onsager Principle 



jump has occurred, the probability that no jump occurs 
in the time interval (0, t) is 



%l){t')dt'. 



(4) 



The waiting time distribution function and the mem- 
ory kernel can be related to one another by taking the 
Laplace transform of the GME ((TJ and the CTRW ©. 
This comparison, after some algebra yields 



itp{u) 
1 - {i^iu) 



(M - I), 



(5) 



where I is the my. m unit matrix and the Laplace trans- 
form of the function / [t) is / {u). Here, as in Allegrini 
et al. we limit our discussion to the two-state case 
where 



M 



1 

1 



and 



(6) 



We approach the subject of fractional derivatives from 
a perspective similar to that of Sokolov and Metzler ■ 
More specifically, we address the problem of the connec- 
tion between the GME |0| and the stationary version of 
the CTRW The GME considered by Allegrini et al. 
(igf is the two-site version of the following equation 



dt 



$(t - t')^^{t')dt' , 



(1) 



where p {t) is the m-dimensional population vector of m 
sites, K is a transition matrix between the sites and ^{t) 
is the memory kernel. The CTRW prescription for this 
process yields 



p (t) = V / dt'i^n it') ^{t- t') M"p 



(0). (2) 



Note that V'n {t) is the probability that n jumps occurred 
and that the last jump took place at time t — t' , implying 
the renewal theory relation 



^nit)= / ^|M-l{t^t')iJl{t')dt', (3) 
Jo 

where ^pi (t) is the waiting time distribution function 
-0 (t) introduced into CTRW and (t) = 5 (t) . While 
M is the transition matrix connecting the sites after one 



K 



1 -1 
-1 1 



(7) 



thereby reducting ^ to 



uip (u) 
l-tp{u)' 



(8) 



This relation between the Laplace transform of the mem- 
ory kernel and the Laplace transform of the waiting time 
distribution function was first obtained by Kenkre et al. 
and is reviewed by MontroU and West [20l |. 
In the case when the lattice has only two sites, 
a left and a right site, the random walker corre- 
sponds to a dichotomous signal ^, with the values 
^{t) = —1, for the left site, and S^{t) — 1, for the 
right site. For the sake of simplicity, we assume the 
two states to have the same statistical weight. Also 
in the two-state CTRW, if we adopt a discrete time 
representation, the motion of the random walker cor- 
responds to a symbolic sequence {^}, with the form 

+ + + + h-t--t--l ....}, 

which shows a significant persistence of both states. 
The waiting time distribution ip (t) is the distribution 
of the patches filled with either -|-'s or -'s. We assume 
symmetry between the two states and a finite first 
moment of ip (t) making it possible for us to define the 
autocorrelation function for the fluctuations ^ {t) 



< mm > 
<e> ' 



(9) 
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because the process is stationary in time |2l'. 

In his original work Onsager considered the case of 
a macroscopic fluctuation that regresses to equiUbrium 
through the phenomenological equations of motion. Here 
we adopt an extension of the Onsager Principle made by 
Allegrini et al. to the case of two states, using the 
probability of the random walker being in state i = 1,2, 
Pi {t), at time t, 



Pi jt) -P2 (t) 

PI {0)~p2 (0)- 



(10) 



The GME cannot provide information about the auto- 
correlation function of the fluctuations, but only about 
the probability of occupying a given state i — 1,2, pi {t), 
at time t. Thus, given that the state is initially out of 
equilibrium pi (0) — p2 (0) ^ 0, the generalized Onsager 
Principle given by (|10|) is the only way to relate the bath 
auto-correlation function to the dynamics of the system. 
Assuming a regression to equilibrium in such a way as 
the retain (|10|l we obtain from the GME ^ using the 
coupling matrix Q: 



dt 



-2 J ^ {t - t') {t') dt' . (11) 



Thus, the Laplace transform of the auto-correlation func- 
tion can be related to the Laplace transform of the mem- 
ory kernel by 



$5 (u) 



2$(u) 



(12) 



The authors of Ref. studied the problem of es- 
tablishing a complete correspondence between the GME 
and renewal theory, characterized by the waiting time 
distribution 



m = (m - 1) 



(T + ty 



(13) 



with /X > 2 to flt the stationary condition. These authors 
determined that this problem could be solved by express- 
ing the CTRW in stationary form, resulting in the GME 
memory kernel: 



(l - ^ {u) 



$ (") = 7 (14) 

-2 [l - ijj {u)j +u[l + iP (m) j r 
where r is the average waiting time 



t4,{t)dt^^. (15) 
fi — z 



The form of the memory kernel given by (|14(l is consis- 
tent with the equation of motion for the auto-correlation 
function l(TT|l . and consequently, Eq. l(T^ is equivalent 
to 



2 $^(u) 



(16) 



We have to remind the reader that the stationary auto- 
correlation function of £ is not related directly to ip{t). 
Zumofen and Klafter |22l | provided a prescription for 
deriving the corresponding equilibrium auto-correlation 
function of £ from ipit)- Their result rests on the observa- 
tion that V'(0 is experimental function, evaluated by 
observing the time duration of the two states. The con- 
nection with renewal theory is established by assuming 
that the time duration of a state is determined by two 
processes, one is the extraction of a random number from 
a theoretical inverse power law distribution, ^*{t), with 
the same power index /x, and the other is a coin toss- 
ing procedure that determines the sign of this laminar 
region. Thus, a given experimental sojourn time in one 
of the two states may correspond to an arbitrarily large 
number of drawings and coin tossings. Renewal theory is 
used to relate the auto-correlation function <i>^(t) to the 
waiting time distribution function V'* {i)- lu fact, from 
the renewal theory |2lj | we obtain the following important 
result: 



1 



it' - t)%l;* {t')dt' , 



(17) 



where t* is the mean waiting time of the '<p*{t)~ 
distribution density. It is interesting to notice that this 
equation implies that the second derivative of the auto- 
correlation function is proportional to ip* (t) , 



(18) 



In Section 2 the departure point of our calculations is 
given by the auto-correlation function <I>^ (t) of Eq. H17|l . 
In this case is convenient to assign to this equilibrium 
auto-correlation function a simple analytical form. This 
is done as follows. We assign the form of Eq. H13|) to 
ip*(t), thereby writing 



{t + T*Y 

This makes it possible for us to write r* as follows 

tr {t)dt = ^. 

p z 



(19) 



(20) 



With the choice of Eq. (fT^ for tp*{t), the autocor- 
relation function ^^(t) of Eq. (|17|) gets the attractive 
analytical form 



$5(i) 



T 



*0 



(T* +t)l3' 



(21) 
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where 



f3 = fi-2. 



(22) 



Thus, in the case < 3, the auto-correlation function of 
the fluctuations is not integrable. 

Zumofen and Klafter |22|, in addition to explaining 
with clear physical arguments the connection between 
ip{t) and '4'*{t), established that the Laplace transforms 
of the two functions are related one to the other by 



first exit times. On the other hand, we extend the ap- 
proach to systems of any age and reveal the phenomenon 
of a continuous time random walk with rejuvenation. To 
accomplish this dual role we rely heavily on the results 
recently obtained by Barkai [^_^and, to a lesser extent, 
the results of AUegrini et al. (13 . However, this allows 
us to reveal some aspects of aging and rejuvenation de- 
pendent on the order of fractional derivatives, which were 
not previously identified. 



riu) 



1 + Tpiu) ' 



(23) 



This important relation allows us to establish a connec- 
tion between r and t*, which turns out to be 



2r* 



(24) 



The equivalence between Eq. (|14|l and Eq. H16|l rests 
on the key property of Eq. H23|l . We note that if we make 
the choice of Eq. (|13|l then the waiting time distribu- 
tion ip* loses the simple analytical form of Eq. H19|l . and 
viceversa. On the same token, the choice of the analyti- 
cal form of Eq. (|13|) for ip{t) makes the auto-correlation 
function $^ (t) lose the analytical form of Eq. H21|) . How- 
ever, using the property of Eq. (|23|l . it is straigthfor- 
ward to prove that ip{t) with the form of (|13|l yields, for 
the auto-ocorrelation function ${(i), the following time 
asymptotic behavior: 



1 



(25) 



Thus, whatever choice is made, either the analytical form 
of Eq. (|13|l or the analytical form of Eq. (|19|l . in both 
cases the two waiting time distributions maintain the 
same time asymptotic behavior, with the same fi. So 
do the two different expressions for the equilibrium auto- 
autocorrelation functions, the time asymptotic equiva- 
lence being the property that matters to study the emer- 
gence of fractional derivative. 

Herein, using the inverse Laplace transform of we 
determine the unknown memory kernel $(t), making it 
possible to discuss how to express the GME in terms 
of fractional derivatives. The case where 2 < /i < 3 
is compared to the recent work of Sokolov and Metzler 
[T^. We find that the index of the fractional deriva- 
tive is 3 — /I, rather than /i — 2, as predicted by Sokolov 
and Metzler. We prove that this difference in index is 
due to the fact that we adopt a stationary condition, 
while Sokolov and Metzler do not. We also prove that 
in the case of a finite, rather than infinite age, our GME 
makes a transition from the (3 — /.t)-th to the (/x — 2)-th 
order. The stationary case becomes stable only in the 
limiting case of infinite age. Thus, on the one hand we 
shed light on the meaning of the work of AUegrini et al. 
|l9j| . which is proven to be a subordination to a Markov 
master equation through the stationary distribution of 



II. THE INVERSE LAPLACE TRANSFORM OF 
THE MEMORY KERNEL 

To establish the form of the unknown memory kernel 
^{t), we make a few preliminary observations. First of 
all, we note that through Eq. H16|) we establish a direct 
connection with the correlation function ^^{t) and that 
this correlation function is, in turn, directly related to 
the waiting time distribution ilj*{t), through Eq. I|17|l . 
Thus, with no loss of generality for the reasons illus- 
trated in Section f , it is convenient to refer ourselves to 
ijj*{t) rather than to V-'(^)- For simplicity, we set T* = 1 
throughout the section. Thus, the Laplace transform of 
the autocorrelation function is Hjl : 



$5 (u) 



r(i-/3) 



/,l-/3 



13-1) 



(26) 



where < (3 < 1, given the fact that we are consider- 
ing 2 < /i < 3, and EJ!j_-^ is the generalized exponential 

function Thus, $5 (u) diverges as u — > and Eq. 
((T^ yields $ (0) = 0. We explore the opposite limit, 
M — > 00 using Eq. which yields ^ (u) ~ ^ — 

In the time representation, the latter limit is equivalent 
to $ it) K, 1^ for t — > 0. Therefore, we segment the 
Laplace transform of the GME memory kernel into two 
parts as follows 



(27) 



The first term models the short-time limit, while the sec- 
ond term is responsible for the long-time behavior. In 
the time representation we have 



It* 



(28) 



Note that this division of the memory kernel into a white- 
noise contribution and a slow term corresponds to a sim- 
ilar partition made by Fulihski [2fil |. 

Thus, for the time evolution equation of the correlation 
function of S,(t), we derive the following equation 



^ _ 1$^ (i) _ 2 1 it - 1') (f) dt'. 



(29) 
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Using Eq. and substituting into it the explicit ex- 
pression of T* as a function of /i, after some algebra, we 
obtain 



T 



t + T J T{t + T) 



it) 



T{t + T) 



(30) 



= <^>a{t~t')^i{t')dt'. 



The two terms on the left-hand side Eq. Ij^UI) are positive. 
Due to the negative sign on the right-hand term of this 
equation we conclude that it might well be that (t) is 
always negative. 

Let us concentrate on the case /3 < 1: using the auto- 
correlation function <i>^ (t) of Eq. (|^ (with T* = 1) and 
using the change of time variable t' + I — > i', we rewrite 
Eq. (1201) in the form: 



t+i 1 

2 J (< + 1 - t') ^dt' + 2 j ^a{t+l-t') ^dt' 



(3t 



13+1- 



(31) 



In the limiting case i — > oo we neglect the second term 
on the left-hand side of this equation. This is a natural 
consequence of the assumption that the memory kernel 
must tend to zero with a negative tail, as an inverse power 
law. With this assumption it is straightforward to prove 
that the modulus of the first term becomes much larger 
than that of the second-term on left-hand side of this 
equation. The consequences of this crucial assumption 
are supported by the numerical results depicted in Fig. 
(^. With this assumption Eq. (|31|l simplifies to 



-2 j^a{t- t') ^dt' = 



(^+1) 



P+1 ' 



(32) 



which can be solved by means of the fractional calculus 
We use the Riemann-Liouville (RL) definition of the 
fractional integral 



1 f f{t')dt' 



(33) 



which is the anti-derivative of the fractional derivative 
with order g, with q < 1. Consequently for /3 < 1 we can 
express (|32|l in terms of the RL fractional integral 



(3t 



2r(l-/5) (t + i) 



0+1 ■ 



FIG. 1; The slow component of the memory kernel $(t), 
"I>a(t), as a function of time. The black dots denote the re- 
sult of the numerical inversion of the expression in Laplace 
transform resulting from Eqs. 11411 and (I27II for /3 = 0.5, the 
continuous line is the analytical approximation given by Eq. 



SO that inverting this equation we have the formal ex- 
pression for the slow part of the memory kernel 



{t) 



2F (1 - P) 



2F (1 - P) 



D]'^ m . (34) 



Pt 



{t + l) 



0+1 



We denote the correlation function with the inverse 
power [P + 1) by the symbol $5,/3+i {t). Carrying out 
the required calculations on the right-hand side of H34|l . 
we obtain, see for example pg.90 of West et al. 2], 



<f a {t) - - 



P 



2T{l~P)T{l + P) (t + iy 



P < 1 (35) 



A comparison with a numerical inversion of the kernel is 
shown in Fig. Q. 

For the sake of completeness, it is worth noticing that 
we can proceed in a similar way also in the case P > 1. In 
this case we find, for the contribution $a(i) of the GME, 
the following time asymptotic behavior 



PiP-l) 



t 



(t+l) 



(36) 



III. 



THE EMERGENCE OF FRACTIONAL 
OPERATORS 



In this section we show that in the two-site case we are 
discussing, the GME has the form of a transport equa- 
tion, with two terms on the right-hand side. The first 
has the form afforded by the ordinary master equation 
and consequently satisfies the Onsager Principle, giv- 
ing a relaxation dependent on the average waiting time 
of CTRW. The second term corresponds to a fractional 
derivative in time, and extends the Onsager Principle to 
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the case of a relaxation with a fat tail. To obtain these 
results, we use what we have learned in the preceding 
section. 

First of all, since we are dealing with the two-site case, 
using the form of M and K matrices, with z = x = —I 
and y = 1, we rewrite Eq. in the following form: 



Ap,(t)=J^it-T){p2iT)-p,iT))dT , (37) 


ip2{t)^J<fit~T){p,{T)-p,{T))dT . (38) 


The memory kernel $ is related to the autocorrelation 
function of the dichotomous variable <I>^ through Eq. 
(|16|l . Inserting Eq. (|16() into the Laplace transform of 
the set of the two-site dynamical equations, solving the 
resulting set of equations, and taking the corresponding 
inverse Laplace transforms yields the solution 

Pi(t) = i[l-$a<)(P2(0)-pi(0))] (39) 
P2(t) = i[l + *?W(P2(0)"Pi(0))] • (40) 

Note that these solutions can be combined to yield the 
generalized Onsager Principle given by H10() in terms of 
the difference in the probabilities. 

We now want to find a formal equation of evolution 
for the probabilities involving fractional operators. We 
know from the preceeding section that: 



$(i) = 



2t* 



*a (t) 



(41) 



with T* — — i thanks to the fact that we set — 
1. Substituting the decomposition of the memory kernel 
into Eq. H37|) . we obtain 



Also we found that in the case < /3 < 1 the asymptotic 
behavior of the memory kernel is expressed by 



$a it) 
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2r(i - /3)r(n- /?) ^t + if 



2r(i - /3)r(i + /?) 

Rewriting Eq. (gSJ as: 

1 



(45) 

t^oo, T* = 1. 



(t) 



2T{1- I3)T{I3){I3 -l)dt 



/3-1 



(46) 



we identify f{t) with 2(i-;3)r(\-/3)r(ff) ^^ ^ for i ^ oo. 
Choosing /(O) — and using the properties of the gamma 
function, we assign to the time asymptotic equation of 
motion the following form 



dpi{t) pi{t)~p2{t) 



dt 



2t* 



(47) 



2r(2 - I3)V{I3) dt 



{t-Tf''{p2{T)-p^{T))dT 



or, in terms of a RL fractional integral l|33|) . 



dpi{t) p^{t)-p2{t) 



Dl"'[p,{t)~p2{t)]. 

(48) 

The same procedure applied to the equation of motion 
for P2{t) yields. 



dt 



2t* 2T{2 - (3) 



dp2it) _ P2{t)~pi{t) 



1 



dt 



2t* 2r(2 - P) 



Dr''[Pi{t)~P2it)]. 

(49) 



The difference between these two equations yields 



dpiit) _P2{t)-pi{t) 



dt 



2t* 



+ J $a(i-T)(p2(r) -pi(T))dT. 


(42) 

Writing (f>a (t) as the derivative of an as yet unspecified 
function /(t), and using the property 



- J f{t-T)g{T)dT^f{Q)g{t) + j f'{t^T)g{T)dT (43) 



with f'{t) = A/(t), we obtain 



dpi{t) _ P2{t)-pi{t) 



dt 



2t* 



f{0){p2{t)-pi{t)) (44) 



+ ^ / /(<-r)(p2(r)-pi(r))dr. 



^^Mt)-pAt)] = -Ei(!1_eM 



(50) 



1 



Dl-^[p,{t)~P2{t)] 



r(2 - /3) 



clearly showing the two kinds of contributions to the gen- 
eralized Onsager Principle. The first term gives the re- 
laxation of the perturbation away from equilibrium at 
the macroscopic rate required by Onsager. The second 
term gives the additional slow relaxation in the form of 
the fractional integral. 



IV. AGING ORDER 

We have to remark that the condition of Eq. p2|) 
refers to the stationary condition explicitly considered by 
Klafter and Zumofen |2^. In this section we prove that 
there is a connection between a system's age and the or- 
der of the fractional derivative in the relaxation process. 
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A sign of the dependence of the fractional derivative or- 
der on age is given by the discrepancy between the result 
of Section 3 and Ref . • Let us compare Eq. (|48|l to 
Eq.(16) of Ref. We obtain the fractional index 1 — P 
rather than (3 as in the work of Sokolov and Metzler. In 
Appendix 1 we prove that our time asymptotic approach 
to fractional derivatives, in the non-stationary case stud- 
ied by Sokolov and Metzler yields the same index as they 
obtain [l^l- Thus, the discrepancy between our predic- 
tion and the prediction of Sokolov and Metzler depends 
on the fact that we consider a condition consistent with 
the Onsager principle, whereas Sokolov and Metzler do 
not. Furthermore, if the system is not infinitely aged, a 
sort of rejuvination process is expected to take place that 
will lead to the fractional order of Sokolov and Metzler. 

To support our remarks concerning the relation be- 
tween aging and the order of fractional operator, here we 
discuss how to define a waiting time distribution of any 
age. The authors of Ref. have shown that the wait- 
ing time distribution tp{t) of Eq. H13|l is obtained from 
the following dynamic model. A particle moves in an the 
interval / = [0, 1] driven by the equation of motion 




FIG. 2: The effective power index as a function of 

time, for /i — 2.3. The curves refer, from top to bottom, to 
ta = 100, 1000, 10000. 



age. Thus, we have that 



lo i^jt + y)dy 

9{ta) ' 



(52) 



where g{ta) is the normalization factor defined by 



dy/dt = ay^ 



(51) 



with z > 1. When the particle reaches the border y = 1, 
it is injected back to an initial condition between y = 
and y = 1 with uniform probability. The age of the 
CTRW is determined by the distribution of first exit 
times. The ordinary CTRW rests on identifying this dis- 
tribution with This means that the CTRW is equiv- 
alent to assuming that the system is prepared in a flat 
distribution at i = 0, which coincides with the beginning 
of the observation process. 

Let us discuss now the consequence of beginning the 
observation a significant time after the preparation. Let 
us imagine that the system is prepared in a flat distri- 
bution at a time t — —ta < 0, and that the observation 
begins at i = 0. This means that the flat distribution 
begins producing a sequence of time intervals of size r, 
according to the distribution of Eq. (|13() ; more precisely, 
the time interval ri beginning at t = —ta and ending at 
t = — ta+Ti, the time interval T2 beginning at t = —ta+Ti 
and ending aX t = —ta + ti + T2, and so on. The waiting 
time distribution of age ta, denoted by 'ipt^{t), is deter- 
mined by the first of these time intervals overlapping with 
t > 0. The time length of that overlap is the time length 
whose distribution determines ^pt^ (t) . We make the as- 
sumption that the beginning of the first time interval 
overlapping with t > occurs with equal probability at 
any point between t — —ta and t — 0. The validity of this 
assumption is discussed in Appendix 2, which establishes 
that this assumption is very good for ta ^ and ta ^ oo. 
In between the asymptotic limits the resulting prediction 
is not exact. However, since it yields simple analytical 
formulas, we adopt this simplifying assumption for any 



gita) = / ^it')dt' 
Jo 



(53) 



and ^'(i) is the probability that no event occurs through- 
out the time interval of length t. Using for ilj{t) the ex- 
plicit form of Eq. ((T^ . it is easy to prove that Eq. l|5^ 
can be written under the form 

/ m ( ,, (t + r)(i-^)~(t + r + tj(i^^) 

This formula proves that for t << ta the index of the dis- 
tribution is /i — 1, whereas for t » ta the index becomes 
/i. This result for the age-dependent waiting time distri- 
bution function agrees with the predictions by Barkai 
and by the authors of Ref. . Notice that the formula 
Eq. I|54|) is equivalent to drawing the initial condition for 
y from an aged distribution of this variable. 

Here, we are in a position to evaluate the waiting time 
index at a generic time, where, we write iptait), as 



V'ta(<) = 



AiT,ta) 

{t + TY-ff 



(*)■ 



(55) 



Using Eq. (|54|l we arrive at the following formula for the 
time dependence of the effective power-law index 

H(t + r)a^M)_(t + T + ^,)(wo]] 

Fig. Q illustrates the regression of the effective power- 
law index to /i with two different ages, and shows clearly 
that the regression is slower for older systems. This for- 
mula does more than explain the discrepancy between 
Eq. igHl and Eq. (16) of Ref . fT5|. In fact, it shows that 
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it is possible to build a GME that at short times follows 
the prescription of our GME and at long times moves 
into the basin of attraction of Sokolov and Metzler. This 
is certainly the case, if ta > —oo. 

This aspect is important and needs a more exhaus- 
tive illustration. We note that the approach of Ref. 
can be easily extended to the case where the distribu- 
tion of first exit times has a finite age. It is enough 
to follow the procedure of Ref. |0 and to replace the 
first exit time distribution with ip^^ (t ) ra ther than with 
ipoo{t) = '0«a=-oo(i)) as done in Ref. jl^l- The result of 
this procedure yields for the GME the following form for 
the Laplace transform of the memory kernel 



(57) 



It straightforward to prove that for ta — Q H57|) reduces 
to 



(58) 



In fact, the general prescription of Eq. 152|l immediately 
yields V'ta(^) — i'i't) foi' = 0. To derive the memory 
kernel corresponding to the infinitely aged condition of 
Eq. H14|) we have to notice first that Eq. (|52|l yields, in 
accordance with Refs. 231 and Il9ll. 



1 



dt'ij{t'). (59) 



To illustrate the change of the memory kernel with 
time, notice that the Laplace transform of Eq. H54|l is: 

. *(M)(l-e"*")-he"*'' f*"e-"2'^'(y)d?/ ^ ^ 



that is: 



9{ta) 



ugita) 



(61) 



Let us substitute H61|l into Eq. H57() . from which we 
obtain the Laplace transform of the memory kernel for 
an arbitrary age. In Fig. Q we show the Laplace trans- 
form of the memory kernel corresponding to a number of 
different ages. Moving from the top to the bottom curve 
the age increases from the brand new kernel {ta = 0) to 
the infinitely aged, or stationary kernel {ta = oo). 



V. CONCLUDING REMARKS 

In this paper we establish that an infinitely old sys- 
tem, with the power-law index in the interval 2 < /i < 3, 
yields the fractional order 1 — /3 = 3 — /i. The prediction 
of Ref. ^2 1 O'^ the other hand, yields the fractional order 



0.6 




FIG. 3: The Laplace transform of the ta-old memory kernel 
^ta{t), ^t„(u) as a function of u, for /i — 2.3. The curves 
refer, from top to bottom, to ta = 0, 0.1, 1, 10, oo 



fj. — '2, corresponding to the brand new condition. If the 
system is not infinitely aged, namely ta < oo, the short- 
time behavior of the system, for < << ta is expected to 
be that of an aged system. At large observation times, 
t » ta, rejuvenation begins. This can be explained us- 
ing the dynamical model of Eq. H51|l . In fact, aging has 
to do with the slow regression to equilibrium, if it ex- 
ists, of the variable y, which mimics a real bath slowly 
regressing to equilibrium. The aging effects discussed by 
Barkai [23 corresponds to z > 2 and thus to /i < 2. In 
this case the dynamic model under discussion does not 
have an invariant distribution, and, consequently, any 
observation is done while the bath is drifting towards a 
condition that will be never reached. This aging effect 
affects the form of the first exit time distribution, whose 
index is — 1 rather than fi. However, after the first exit 
the trajectories are injected back with a uniform proba- 
bility, and thus, all the ensuing jumps are determined by 
the ordinary waiting time distribution 'ijj{t). 

It would be desirable to have an equation of motion 
with a fractional operator order that changes as a func- 
tion of time from 3 — /i to /i — 2. 

However, there are technical and conceptual difficul- 
ties that make it difficult, if not impossible, to realize 
this goal. In fact, according to the perspective adopted 
in this paper, the order of the fractional operator is es- 
tablished using time asymptotic arguments. Thus, if 
< oo we can associate a time t « ta with the or- 
der ord = /i — 2, of the fractional derivative. This is so 
because for t « ta the calculations would be virtually 
equivalent to those done in Appendix I. It is not clear how 
to proceed when t is is of the order of ta, this being the 
first reason why assigning a fractional derivative order to 
any time might be difficult. There also exists a physical 
reason that might make it impossible to move from the 
order 3 — /i to /i — 2. Physically this extremely extended 
transition process might involve a mixture of fractional 
derivatives of different orders. It has been shown |23| 
that the Levy walk does not have well-defined scaling, 
due to aging effects. Similarly, the adoption of a frac- 
tional time derivative, with time-dependent order, might 
be inadequate to explore the regime of transition from 
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the order 3 — to the order ^ — 2. In conclusion, the 
fractional operator, and its order, reflects a stable condi- 
tion, of a brand new or infinitely old system. The regime 
of transition from the dynamic to any of these two ther- 
modynamic regimes, and the regime of transition from 
the earlier to the latter thermodynamic regime is not yet 
a fully understood physical condition, an issue calling for 
further investigation. 

It is interesting to notice that, even if we select n > 2, 
and consequently we adopt a condition compatible with 
the stationary condition, the effective index of the first 
exit distribution is located in the non-stationary region, if 
/3 < 1. This is probably the reason why the memory ker- 
nel seems to share the same properties as those adopted 
by Ref.[^ ITol l28| to produce sub-diffusion. Notice that 
the baths used by Lutz and Pottier have proper- 
ties quit e different from the subordination perspective of 
Ref. ha , even the relaxation process stemming from sub- 
ordination ^3 is quite similar to that produced by the 
non-ohmic baths of Lutz and Pottier. We hope that the 
present work might help understanding the connection 
between the two perspectives. This is another subject 
for future research. 

Acknowledgments G.A and P. G. gratefully ac- 
knowledge financial support from ARO through Grant 
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stituting it into Eq. ljA-l|) we have: 

$ (M-imi-M)(e"-j^;_,) 
^x{^i - l)T (-/i) u (e" - El) 

u^-^^-{^Ji-l)T{l~^Ji) (e--i?-_i)' 

where, as usual, for simplicity we have set T — \. An- 
ticipating the convolution form of the solution we cross- 
multiply to obtain: 

[u'-^^ - \)T{\-^i) (e" - El_^)] {u) 
= (/.-l)((M-l)r(l-M) 

(A-3) 

Using the following relation 0] 

oo 

/ exp \-ut] dt = iEl'' - e""^) , a > -1 



(A-4) 

and setting a — fj, ~ 1 we construct 



APPENDIX I 



(A-5) 



In this Appendix we show how to obtain the order 
of the fractional operator in the GME for 1 < < 3, 
using the Laplace transform form for $(<;) given by Eq. ( 
IH)). We are referring ourselves to ip{t) rather than 'ip*{t). 
Thus, for the reasons pointed out in Section 2, we adopt 
the analytical expression of Eq. (|13|) . 

We note that for it ^ oo we get a finite value: 
4 (oo) — corresponding to ip (0) in t-space. As done 
in Section 2, we separate the kernel into two contribu- 
tions: ^{t) = ii^S{t) + ^ait), and insert the Laplace 
transform of separation into Eq. to write: 



(A-1) 



uj^ju) - (/i - 1) + (/i - 1)^(m) 
1 - tfjiu) 

iJ'D ju) + (/i - l)-0 ju) 

1 - (u) 



We introduce 4'd (u) is the Laplace transform of the dis- 
tribution's derivative ip' (t) = —fJ-ifJ- — 1) (^j^iy+i ; using 
the Laplace transform of an inverse power law and sub- 



where R{t) is the inverse Laplace transform of right side 
of Eq. (|A-3p . Using the well known recurrence relation 
of Gamma function, we can combine terms in Eq. (|A-5|I 
to obtain 



m - / 1'"-' - ) Mt - t')dt' 







(A- 



Let us consider first the case > 2, where, in the limit 
of i ^ 1 (equivalent to t ^ T), we can write Eq. (|A-6|) 
to a good approximation 



r 

R{t) = fr^j^—Yj J t'^^'^ait - t')dt' = R{t). (A-7) 
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Going back to the Laplace transform representation, we 
obtain the simplier expression 



r(^-i) =(A^-i) r(i-M)(e 

-/.0.-i)r(-M)^/(e"-i?;'), 



(A-8) 



which after a Httle algebra yields for $a (u): 



$a (m) = - 1) - 2) 



(M-i)r(i-Ai)- 



1 e - 



u u 



1 e" - 



(A-9) 



So that using the inverse Laplace transforms we obtain 
the corresponding expression in the time representation: 



{t) - (a* - 1) (a* - 2) 



In the time asymptotic limit we get 



(A-10) 



overlapping with i > is uniformely distributed between 
t = —t„ and t = Q. Let us set t„ = 0. 



The exact expression for ij^t^ (t) is 



il^tM = / dxG{ta-x)^{t + x), 
Jo 



(B-1) 



where 



Git) = S{t) + ij{t) + V / driV(Ti) (B-2) 

t /.t 
dT2lp{T2 - Tl)... dr„_l7/)(t - T„_l). 




1 



(A-11) 



corresponding to fractional operator of index f3 = fi — 2. 
For the sake of completeness, we also give the expression 
for the GME kernel in the case /x < 2. Proceeding as 
done earlier we obtain: 



FIG. 4: The waiting time distribution ij^tait) as a function ot 
time t. The dots denote the exact values, and the lines the pre- 
diction of Eq. I|54|l for ^ = 2.3, T = 1. Moving from bottom 
to top in the right-hand portion of the ta ~ 0.01, 1, 10, 100 



sinTT/i t'' 



[^l-^ + i^i-2)t]. (A-12) 



Finally, we want to point out that the expressions we are 
proposing refer to values of mu which are not integer. We 
are exploring the interval [1,3]. Thus the the expressions 
we are proposing become questionable for = 2. To get 
the proper expression for fi = 2 we have to apply to study 
expressions like those of Eq. ljA-l|l and (|A-8|) at /i = 2-f e, 
do a Talylor series expansion around /i = 2, and assign 
to fi = 2 the limiting values reached for e ^ 0. 



APPENDIX II 



This Appendix is devoted to establishing the accuracy 
of Eq. (|52|) . and consequently the validity of the as- 
sumption that the beginning of the first laminar region 



It is straightforward to find the Laplace transform of 
G{t). This is given by 



G{u) = J2i>{ur^- : 



(B-3) 



Thus, the Laplace transform of (jB-l|) respect to ta reads: 



A At) 



1 



l-i'iSa) 



iJKSa) 



-Say 



ip{v)dy 



(B-4) 



By the numerical anti-Laplace transforming Eq. (|B-4p we 
evaluate the time dependence of the exact waiting time 
distribution of age ta, Eq. (|B-1|I . In Fig. I^J we com- 
pare the exact prediction, evaluated numerically, to the 
heuristic expression of Eq. 1)52(1 . We find that at short 
and large values of ta these two expressions coincide. In 
the intermediate region they do not. Nevertheless, we 
think the agreement between the two expressions is sat- 
isfactory enough for the purpose of this paper. 
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